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Q>^ ■ Abstract 

A QCD based effective action is constructed to describe the dynamics of confinement 
and symmetry breaking in the process of parton-hadron conversion. The deconfined quark 



and gluon degrees of freedom of the perturbative QCD vacuum are coupled to color singlet 
collective fields representing the non-perturbative vacuum with broken scale and chiral 



■ symmetry. The effective action recovers QCD with its scale and chiral symmetry properties 

at short space-time distances, but yields at large distances (r >, 1 fm) to the formation of 
symmetry breaking gluon and quark condensates. The approach is applied to the evolution 
of a fragmenting qq pair with its generated gluon distribution, starting from a large hard 
scale Q^. The modification of the gluon distribution arising from the coupling to the non- 
perturbative collective field results eventually in a complete condensation of gluons. Color 
flux tube configurations of the gluons in between the qq pair are obtained as solutions of 
the equations of motion. With reasonable parameter choice, the associated energy per unit 
length (string tension) comes out ~ 1 GeV/fm, consistent with common estimates. 
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1. INTRODUCTION 



The physics of QCD exhibits different relevant excitations at different length scales. At 
space-time short distances (below 1 fm) the relevant degrees of freedom are quarks and 
gluons whose interactions are well described by perturbative QCD Q. The long distance 
physics on the other hand is governed by the hadronic degrees of freedom, and the particles 
which are observed at large scales are hadrons whose interactions are well described by 
chiral models j^. The change of resolution of our microscope with which we probe the 
physics of QCD is formally described by a renormalization group equation, or evolution 
equation, that determines the scale dependence of the theory |@, |5|, |6|. 

The transition from short distance (high momentum transfer) regime to the long dis- 
tance (low energy) domain can be cast in terms of an evolution equation for an effective 
QCD action that embodies both, fundamental partonic degrees and hadronic degrees of 
freedom. By increasing the distance scale (decreasing the momentum scale), the evolution 
1^ of the effective field theory must lead from one set to the other set of degrees of freedom. 

Experiments on high energy QCD processes, such as e~^e~ annihilation, deep inelastic ep 
scattering, Drell Yan, etc., strongly support the conception that the observed parton frag- 
mentation into hadrons is a universal mechanism. Moreover, the dynamical transformation 
of color charged quarks and gluons in high energy QCD processes into colorless hadrons is 
commonly believed to be a local phenomenon ||8|. Thus, a consistent description of the lo- 
cal hadronization mechanism must be independent of the details of the partons prehistory 
and should in principle apply also to hadron-hadron, hadron- nucleus, or nucleus-nucleus 
collisions. 

To date most of the theoretical tools to study properties of QCD are inadequate to 
describe the dynamics of the transformation from partonic to hadronic degrees of freedom: 
Perturbative techniques are limited to the deconfined, short distance regime of high energy 
partons |^, QCD sumrules [10| and effective low energy models |ll[| are restricted to the 
long distance domain of hadrons, and lattice QCD lacks the capability of dynamical 
calculations concerning the quark-gluon to hadron conversion. On the other hand, phe- 
nomenological approaches to parton fragmentation [13|, are mostly based on hadronization 
models with adhoc prescriptions to simulate hadron formation from parton decays. 

In this paper I follow a rather different, universal approach to the dynamic transition 
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between partons and hadrons based on an effective QCD field theory description, as re- 



theory from high energy to low energy scales 0], the key element is to project out the 
relevant degrees of freedom for each kinematic regime and to embody them in an effective 
QCD Lagrangian which recovers QCD with its scale and chiral symmetry properties at high 
momentum transfer, but yields at low energies the formation of symmetry breaking gluon 
and quark condensates including excitations that represent the physical hadrons. In Sec. 
2, I will first formulate the general field theoretical framework. On the basis of the dual 
vacuum picture of coexisting perturbative and non-perturbative domains an effective action 
is constructed that embodies the correct scale and symmetry properties of QCD. The con- 
cept is here more phenomenologically motivated than the related formal approach of Ref. 
[0]. However, there appears to be a clear correspondance bertween those two descriptions. 
In Sec. 3, I shall demonstrate the applicability of this effective QCD field theory to the 
dynamics of parton-hadron conversion by exemplarily considering the evolution of gluons 
produced by a fragmenting quark-antiquark pair. The change of the gluon distribution in 
the presence of a confining composite field is studied and flux tube solutions of the gluon 
field resulting from the equations of motion are analyzed in terms of the string tension, that 
characterizes the effective confinement potential. Various perspectives of the approach are 
discussed in Sec. 4, in particular the applicability to the QCD phase transition, and high 
density QCD. 

2. EFFECTIVE QCD FIELD THEORY WITH SPONTANOUS SYMMETRY 
BREAKING 

The goal is to construct an effective field theory that describes the dynamics of both 
partonic and hadronic degrees of freedom and their interplay. The approach is based on 
the concept of the effective action |^, ^, which will be represented here as (r = denotes 
the space-time 4-vector) 



The three contributions to the action, which will be discussed below, correspond to the 
QCD Lagrangian with the quark (ipi) and gluon fields {Aa), an effective low energy La- 
grangian introducing composite fields x and U, and a term that couples the "microscopic" 



cently proposed in Ref. 1 14 1 . In the spirit of the aforementioned evolution of effective field 




(1) 
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fundamental quark and gluon degrees of freedom to the " macroscopic" fields x ^-iid U which 
represent the hadronic degrees of freedom. 



2.1 Clip, A] 

The QCD Lagrangian in (|l]) contains the gluon fields coupled to massless quark 
fields = 1) • • • ) ^f)i 



^) + ^}Hn}l + ^}(i/u • (2) 



Here F>^'^ = A^^ — + ggfabcA'^A'^ is the gluon field strength tensor. The subscripts 

a,b,c label the color components and denotes the color charge related to as = 5'g/(47r). 
The Ta are the generators of the SU{3) color group, satisfying [T^jTb] = ifabcTc with the 
structure constants fahc- The gauge fixing term >C}^n}] = oo/(€H)(r7^^'^)^ with gauge 
parameters a and rj^, and the contribution of Fadeev-Popov ghost fields C; ^}{i!u = 
{dp,C,*\){5^\^d^ — will be irrelevant later on, because a physical gauge -q ■ A = d 

can be fixed, which eliminates the presence of ghosts. 

The Lagrangian (Q) is well known to be invariant under chiral transformations [Q. 
At the tree level it is also invariant under scale transformations ^ r'^ = e^r^ |15], 
generated by a so-called dilaton charge D{t) = / d^rJ^^r), where is the scale current 
and [ D, ip{r)] = i [r^d^ + d^) <f{r) for a generic quantum field if with scale dimension d^. 
The convention is (i^ = 1 for gauge boson fields and d^ = 3/2 for fermion fields. It follows 
that C{'ip,A) has scale dimension 4 so that [D,C{ip,A)] = and therefore massless QCD 
proves to be scale invariant at tree level. 

At high energies and short space-time distances, asymptotic freedom leads to unconfined 
gluon and quark fields in (^). However, in the physical world these color degrees of freedom 
are confined, and both chiral and scale symmetry are explicitely broken. To describe the 
dynamics of the symmetry breakdown of the transition between the perturbative, scale 
and chiral invariant, regime and the non-perturbative world with broken symmetries, one 
needs to supplement (Q) (by adding AC = C[xM\ + x]) to construct an effective 

description such that at high energies the fully symmetric QCD phase is recovered, but at 
low energies massive hadrons emerge. 
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2.2 C[xM] 

The specific form of C[xM\ in (H) is adopted from Refs. |16, 17], where an effective low 
energy Lagrangian was constructed, guided by the scale and chiral symmetry properties 
of the QCD Lagrangian. The construction is based on the observation that even massless 
QCD is no longer scale invariant when going beyond tree level, because of the scale anomaly 
[p^ d^J^ = f3{gs) / {2gs)F^uF^^ , where is the scale current as before. In addition chiral 
invariance breaks down when finite quark masses are taken into account. As a consequence, 
the QCD energy momentum tensor 6f_iy exhibits the well known trace anomaly [p^ ]: 

/3(a 



6*^ 



4a« 



F.^F^^'^ + (1 + 7^) E 



m. 



(3) 



where (3 is the Callan-Symanzik function and 7^ is the anomalous mass dimension. This 
anomaly constrains the form of the effective low energy Lagrangian, because without it 
Poincare invariance would be broken, and the mass of the proton would come out wrong, 
since 2m'^ = (^16*^1^). 

The extension of these symmetry properties of the QCD Lagrangian to the low energy 
domain was modeled by Campbell, Ellis and Olive [0] asQ 



\{d^x){d^x) + \Tr 



{d^U){d>'U^) 



4 , 4 1 



X 



Xo 



c Tr 



mg{U + U^) 



x_y 

XoJ 



1 2^2 fxY 



(4) 



This form introduces a scalar gluon condensate field x and a pseudoscalar quark condensate 
field U = fiT exp (j,J2^=o ^j4'j / Itt^ for the nonet of the meson fields (pj (/^r = 93 MeV, 
Tr[XiXj] = 26ij, UU^ = f^), with non- vanishing vacuum expectation values 



Xo 



(o| x|o) / , 



c{0\U + \0) / 



(5) 



that explicitly break scale, respectively chiral symmetry. In (^), b is related to the con- 
ventional bag constant B hy B = 6xq/4, c is a constant of mass dimension 3, ruq = 
diag(mu, m^, ms) is the light quark mass matrix, and tuq is an extra U(l)-breaking mass 
term for the ninth pseudoscalar meson (f)Q. 

^In Ref. the second term in eq. was scaled by a factor (x/xo)^, which however is not necessary 
Pq], since just this term gives a chiral symmetric contribution in agreement with the QCD anomaly (0). 
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Notice that the anomaly constraint is modeled by the third and fourth term in (^) 
with the correspondence 

(Ol^F^.F'^nO) = -bxt (6) 



and 



{0\qq\0) = c (—Y {0\U + U^\0) . (7) 
\Xo/ 



2.3 C[^P,A,x] 

The key ingredient in (|l|) is the connection between the scale and chiral symmetric, short 
distance regime of colored fluctuations and the world of colorless hadrons with broken 
symmetries. It is mediated by the coupling between the fundamental quark and gluon 
degress and the collective fields x and U through coupling functions g{x) and ^(x), 

£[V,Ax] = ^ ^MM-^r - V^)}(X)^) ■ (8) 



The coupling functions £,{x) and g{x) are chosen in the spirit of Friedberg and Lee |19], 
who formulated a dual QCD vacuum picture: High momentum, short distance quark-gluon 
fluctuations (the perturbative vacuum) are embedded in a collective background field x 
(the non-perturbative vacuum), in which by definition the low momentum, long range 
fluctuations are absorbed. Confinement is thus associated with the colordielectric stucture 
of the QCD vacuum. This property is modeled by a colordielectric function 

<x) = 1 - e(x) (9) 

that satisfies 

k(0) = 1 , k{xo) = , (10) 

thereby generating color charge confinement, because a color electric charge creates a dis- 
placement Da = nEa, where = F^^, with energy ^Jd^rD^/K which is infinite for 
non-zero total charge if k falls off faster than l/y^ for large distances r. The particular 
form of k{x) is not crucial as long as the properties (10) are satisfied [20|. A specific choice 



is |1| 



k(x) = 1 + ^ 3^-4 e{x), (11) 
\XoJ V Xo 
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which has the further properties of k'(0) = k"(0) = 0. Other forms used in the hterature 
are e.g. k[x) = |1 ~ (x/xo)"!*" [^H (Friedberg and Lee originally proposed n = m = 1). 

Similarly, absolute confinement can be ensured also for quarks by coupling the quark 
fields to the x-field through |23] 



aix) 



90 



1 



1 



(12) 



which leads to an effective confinement potential with the masses of the quarks inside 
approximately equal to the current masses, but at x = Xo it generates an infinite asymptotic 
quark mass (the value of go is irrelevant in the present paper). 

2.4 Equations of motions 

To summarize to this end, the complete effective action (||) is determined by the La- 
grangian 



C 



1{{ 



1 



ilt.d^' - gs7f.A^Ta - g{x) 



l{d^x){d^x) + \Tr 



(5^[/)(a'^C/t) 



v{x,u) , 



(13) 



plus the terms C}-\n}'\ ^"^^ "^Ki/u of (H)- The potential V is given by 



v{x,u) 



1 



h 

+ c Tr 



4X0 + X ln^^,/4 



Xo 



m, 



3 I 

Xo/ ^ 



-)■ 

Xo/ 



(14) 



which has its minimum when x = Xo = (0|x|0) and equals the vacuum pressure (bag 
constant) B = 6xo/4 at x = 0. Typical forms of V{x, U) for different values of B and niq 
are depicted in Fig. 1. 

The effective field theory defined by (|l|) and (^) represents a description of the duality 
of partonic and hadronic degrees of freedom by coupling the high energy QCD phase with 
unconfined gluon and quark degrees of freedom to a low energy QCD phase with confinement 
and broken chiral symmetry which contains a gluon condensate (^ and a quark-antiquark 
condensate Small oscillations about the minimum of the potential V{x,U) are to be 
interpreted as physical hadronic states that emerge after symmetry breaking. They include 
[17|: (i) glueballs and hybrids as quantum fluctuations in the gluon condensate xO) (ii) 



pseudoscalar mesons as excitations of the quark condensate Uq, (iii) the pseudoscalar flavor 
singlet meson <pQ, and (iv) baryons as non-top ological solitons p^ . 

The field equations which derive from (||) and (13) are: 



[7^ (^5^ - gsA^.Ta) - g{x)] i'^ = (15) 



Hx) FH = -9s i^ix) fabc A^^bFr + 9s YT^ A (16) 
d V{x,U) ^ dK{x ) 
dx dx 



df^d^X + ' + -^F^^^aFr + 9ix)i^ii^i = (17) 



dVix,U) 
dU 

Notice that the [/-field does not couple directly to the quark and gluon fields. Per con- 



d.d^U + ' = . (18) 



struction |17], the dynamics of the quark condensate field U is solely driven by the gluon 
condensate field x- It is important to realize that the interplay between the x-field and the 
quark and gluon fields, tp and A, is the crucial element of this approach. 

2.5 Comments 

The following remarks concerning the effective Lagrangian (^3|) are important: 

a) At short distances or high momentum transfers the exact QCD Lagrangian (|2|) is 
recovered, since x = = and k(x) = 1 (i.e. ^(x) = 0) and g{x) = 0, whereas the 
long distance QCD properties emerge as x/xo — ^ 1 and U /Uq — > 1 |l^ and no colored 
quanta survive. The transition from one set of degrees of freedom [ip, A) to the other 
(x, U) corresponds to consecutively integrating out all colored quantum fluctuations 
and absorbing them effectively in the collective color singlet fields. 

b) The problem of double counting degrees of freedom has to be carefully inspected. 
Although it does not arise in one-loop calculations (to which I will restrict here), 
processes with e.g. two-gluon exchange could also be contained in the exchange of 
a color singlet X'Quantum. A minimal possibility to avoid this problem is a rigid 
separation of high and low momentum modes, by introducing a characteristic scale 
Qq: above Qo the physics is described in terms of quark gluon degrees and below Qq 
the dynamics is governed by the collective degrees of freedom Q. 

c) C[xM] for the composite fields embodies the correct QCD scaling and chiral properties 
and accounts for the important anomaly (^ of the physical energy-momentum tensor 



of QCD. The coupling between quarks and gluons to the composite field x in ^[i'-, x] 
can be interpreted in analogy to a thermodynamic system in equilibrium with a heat 
bath, with a net flow of energy between the system and the heat bath environment 
such that the hare energy of the system is not conserved. However the free energy of 
the system, here high momentum quarks and gluons, is constant [^]. It corresponds 
to the conserved energy momentum tensor O^j^y with its non-zero trace (^). 

d) There is no need for explicit renormalization of A£ = C[x-,U\ + £["0,^, x]- The 
composite fields x ^i^cl U are already interpreted as effective degrees of freedom with 
loop corrections implicitely included in A£ and it would be double counting to add 
them again. Moreover, in the present approach the low energy domain of U] is per 
construction bounded from above by the onset of the high energy regime described by 
£[7/^,^4]. In correspondance to item b) the scale Qo that separates the two domains, 
provides an "ultra-violet" cut-off for C[xM], and at the same time an infra-red cut-off 
for C[ip,A]. 

This effective field theory approach offers a wide range of physical applications and can 
be extended and refined in various directions, as discussed in Sec. 4. The scope of the 
remainder of this paper is however concepted as a exemplary demonstration of how of the 
dynamics of parton-hadron conversion emerges within this framework. 

3. CONFINEMENT OF GLUONS IN A FRAGMENTING qq SYSTEM 

The effective QCD field theory defined by (13) is readily applicable to describe the 



dynamic evolution from perturbative to non-perturbative vacuum in high energy processes. 
In accord with the symmetry breaking formalism of Sec. 2, the parton-hadron transition can 
be visualized as the conversion of high momentum colored quanta of the fundamental quark 
and gluon fields into color neutral composite states that are described by the condensate 
fields X s-iid U and their excitations. 

In the following I shall consider as an example the fragmentation of a jet system with 
its emitted bremsstrahlung gluons and describe the evolution of the system as it converts 
from the parton phase to the hadronic phase. The process is illustrated in Fig. 2: A 
time-like virtual photon in an e"'"e~ annihilation event with large invariant mass S> 
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is assumed to produce a qq pair which initiates a cascade of sequential gluon emissions. 
(Here and in the fohowing A denotes the fundamental QCD scale). The early stage is 
characterized by emission of "hot" gluons far off mass shell in the perturbative vacuum. 
Subsequent gluon branchings yield "cooler" gluons with successively smaller virtualities, 
until they are within Qq, where Qo is of the order of = 46xq ~ 1 GeV. At this point 
condensation sets in, or loosely speaking, the "cool" gluons are eaten up by the color neutral 
gluon condensate field Xj the particle excitations of which must then decay into physical 
hadrons by means of some local interaction in the non-perturbative vacuum. A similar 
picture holds for the [/-field. 

It is well known that the bulk of produced particles stems from rather soft gluon emis- 
sions that are characterized by small values x of the fraction of the initial energy. Secondary 
production of qq pairs is comparably rare on the perturbative level. It is therefore reasonable 
to neglect the quark degrees of freedom and study the purely gluonic sector. Furthermore it 
is convenient to work in a physical (axial) gauge for the gluon fields, 77 • ^ = 0, by choosing 
the gauge vector in (§) as r/^ = with n being space-like and constant, in which case the 
ghost contribution vanishes. Consequently, the equations of motion (p^-(p!8D reduce to: 



dA^ix)Fr] = -gs^{x)fahcA^,,tFr (19) 

dVjx) _ di4x) 
dx dx 



d.d^X = - ^F,.,aFr (20) 



X^' 



d^df'U = -4c (^^j Tr[mq\. (21) 

The solution of these equations is a still formidable task, because not only the gluon fields 
but also X ^iid U are quantum operators. To make progress, I will now proceed by (i) 
treating the quantum gluon fields perturbatively, and (ii) employing the mean field approx- 
imation for the composite fields x ^-iid U . Representing 

X{r) = X{r) + X{r) , (22) 

where x is a c- number and x & quantum operator (similarly for U), the mean field approx- 
imation is obtained by neglecting the quantum fluctuations x ™d keeping only x- Thus, 
the approximations (i) and (ii) correpond to the semiclassical limit in which gluonic quan- 
tum fluctuations interact with a classical mean field. To a good approximation this should 
provide a reasonable description: first, because renormalization group improved QCD per- 



turbation theory allows for an accurate description of the evolution of the gluon field [26| at 



short distances where x = 0) &iid second, because the dynamics of the system around x = X 
is governed by a large number of virtual excitations, corresponding to coherent modes of 
field quanta, so that a quasiclassical mean field description should be applicable in the low 
energy regime [19|. 



3.1 Evolution of the gluon distribution in the presence of the collective field x 

As I will show now, the equations of motion (p^)-(21) simplify to a perturbative evolution 
equation for the gluon distribution which is coupled to the equation for the mean field x- 
The key problem is the first equation, since the gluon fields A^, or equivalently F'^'^, drive 
the dynamics of the x-field which in turn feeds back via k{x)- As mentioned before, the U- 
field does not couple directly to the gluon field. The procedure in the following is therefore 
to 'solve' eq. (^) as a function of k and then to insert the solution into (pO|), so that one 
is left (aside from the simple third equation) with a single equation for Xi which however 
is non-linear. 

Solving the equation of motion (|l^) for the gluon field is equivalent to the calculation 
of the complete Greens function with an arbitrary number of gluons. Instead, I will restrict 
to evaluate the 2-point Greens function only (Fig. 3), i.e. the full gluon propagator which 
includes both the one-loop order gluon self-interaction through real and virtual emission 
and absorption, and the effective interaction with the confining background field x- In 
the framework of "jet calculus" [^, this gluon propagator denoted as Dg{x,k'^;xo,Q'^), 
describes how a gluon, produced with an invariant mass Q^, evolves in the variable x 
(momentum or energy fraction) and the virtuality k'^ (or transverse momentum k'j_ ) through 
these interactions. To one-loop order, it is obtained by calculating the corresponding cut 
diagrams. In the present case, one has in addition to the usual gluon branching and 
fusion processes, g ^ gg and gg — > g, contributions from energy transfer and two gluon 
annihilation processes g ^ gx ^^'^ gg ~^ X> respectively 0. This is illustrated in Fig. 3. 

To write down the determining equation for the gluon propagator Dg, it is convenient 
^In the present case of qq jet evolution, the contribution of perturbative 2-gluon fusion processes gg ~* g 
for > is very small [ psj , but the non-perturbative gluon recombination gg ^ x ii^i the range 
is of essential importance in order to achieve complete confinement. 
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employ lightcone variables, defined by the identification of components of four-vectors as 

fc^ = {k+,k-,k^) ; = -^{k^±k^) ■ k^ = {k\k'^) . (23) 

v2 

The k^ component of a particle's momentum, the light cone momentum, is always positive 
definite, k~^ > 0, and the light cone energy, k~ = {k'^ + iv?)/{2k^) is also positive. Further- 
more, the light cone time r+ = (i + z)/\/2 is conjugate to k~ and the light cone coordinate 
r~ is conjugate to k~^ . The invariant momentum space element is 

d^k d^k 2^ dk+d'^k 



Choosing the light cone gauge for the gluon fields, rj ■ A = = —A~ = 0, results in 
well known simplifications in the perturbative analysis of light cone dominanted processes 
and has the advantage that there are neither negative norm gluon states nor ghost states 
present |2^. As a consequence only the transverse components A^j_ (i = 1, 2) are dynamical 
field variables, since A'^ is identically zero and A~ is determined at any 'time' r"*" by 
and . The particular choice r] = [pq + pq) / 2 has the advantage that interference terms 
do not contribute to leading log accuracy (they are suppressed oc 1/k'^). Therefore, in 
the leading log approximation (LLA) [|^, ^, |2^, it is enough to realize that for every 
choice of h (Fig. 3 top), one can group the other gluons in a unique way to groups forming 
dressed rungs of a ladder (Fig. 3 middle) whose discontinuity is taken (Fig. 3 bottom). 

Introducing the variable x = k^/P^ (the light cone fraction), and parametrizing the 
momenta of initial quark and antiquark as P = Pq + Pq with 

Ath?' 

= Q , = ^ ' = ' (25) 

i.e. P^ = Q^, where as before Q denotes the invariant mass of the time-like photon that 
creates the pair, the determining equation for the gluon propagator Dg{x,XQ;k'^,Q'^) can 
now be represented in the form (c.f. Appendix A): 

Dg{x,k^;xo,Q^) = xo6{x-xo)6{k^ -Q^) FgiQ^Ql) 

+ FgiQ^k') [' ^wix',x,e) Dgix',k'';xo,Q') Fgik'^Ql) 

JQl k ^ Jx x' 

(26) 

This equation has a simple physical significance: The first term is the inclusive sum of 
virtual emissions and reabsorptions, and therefore does not change the number of gluons in 
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the gluonic wavefunction of the fragmenting QQ pair, whereas the second term describes 
the change of the gluon distribution as a result of real decay or fusion processes. The 
Sudakov formfactor 



FgiQ^,k^) = exp 



(27) 



is the probability a gluon propagates like a bare, non-interacting particle while degrading 
its virtuality from to k'^. As the gap between and k"^ grows, such a fluctuation 
becomes increasingly unlikely. The total interaction probability 



Wgik^ 



dx I — - w[x ,x,k) 



10 Jx X 

is the integral over the inclusive probability for all possible gluon interaction processes i 



(28) 



proceses i 



w{x',x,e) = J2 . (29) 

The normalization is such that 

rQ^ dk' 2 

1 = Fg{Q\Ql) + FgiQ^k-') / n,g{k;^)Fg{k!\Ql) (30) 

in accord with unitarity (probability) conservation. Multiplying (^) by F^"^, differentiat 



ing, and accounting for (|30|), yields: 

<9 „ , . n „r,, dx' 



k 



d dx 

Dg{x,k'^;xo,Q'^) = / —w{x',x,k'^)Dg{x',k'^;xo,Q'^)-Wg{k'^)Dg{x,k'^;xo,Q'^) 

(31) 

As summarized in the Appendices, one obtains for the individual interaction probabil- 
ities to one loop order (with the assignment xi — > X2, (xi — X2) for branchings and 
xi, {x2 — xi) X2 for fusions) 

Wg^gg{Xl,X2,k'^) 



a. 



. X2 



Wg^gg{Xl,X2,k'^ 



a. 



w 



a-*ax 



'xi,X2,k'^ 



27r 



STT^ Cgg^g Xl{X2-Xl) 

A;2A2 



X2 



Xi 

7™ I - 



27T ^™Ul 



Wgg^^{Xl,X2, k^) 



2tt 



87r2 c, 



xi{x2 - Xl) 



"^x^aa 



(32) 



where c, 



93^9 ~ "^aa^x 



1/8, and 



Ig^ggiz) = 2C A Ul - z) + ^ + 



1 - Z 
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1 ( l + z^ ' 
4 i 1 - z 



z + 



(33) 



Here Ca = Nc = 3, and z is the fraction of x-values of daughter to mother gluons. In 



b In 



A2 



llA^, - 2N 



f 



(34) 



is the one loop order QCD running couphng (in the present case Nj = 0), and 



(35) 



denotes the couphng to the X'^sld in momentum space, with denoting the Fourier 
transform of ^,{x) iii ®- 

Eq. (p6|) for the propagator Dg{x, k'^;xo, Q"^) corresponds the evolution equation for the 
gluon distribution g(x, k\), which is defined ||3^ as average number of gluons at 'light cone 
time' r+ = in the multi-gluon state^ \P) with light cone fractions x = /P^ in a range 
dx and transverse momenta in a range d?k±: 



1 



xg{x,k^) = — 



dr d"^ 



^^^) (P| F+'^(0,r-,f^) F,,^+ (0,0,0^) IPI36) 



(27r)3 

Thus, as is evident from Fig. 3, the probability for finding a gluon with x and k'j_ is given 
by the identification 



g{x,kj 



Dg{x,k^;l,Q'') 



(37) 



On account of the explicit expressions (|32D-(p5|), and taking as evolution variable the gluon 
transverse momentum k'j_ rather than the invariant mass k'^ [^], one finally arrives at the 
master equation for the evolution of the gluon distribution: 



ki 



d 



gix,k1) 



+ 



as{k 



2tt 
2tt 



dz 



-g[-,k_i 

z \z 



99^9 



dz 



.(2) 



7™(^) 0{kl-Ql) 
1 - z 



^5(2) (^,^^ (1 _ kl)] Jg^gg (z) 9 [kl - 



■^It is convenient to visualize the initial qq pair (^5|) as a single incoming 'gluon' with momentum P, i.e. 
with xo = 1 and invariant mass Q (c.f Fig. 3). 
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27r 



dz 







(38) 



27r 



99^x I dz ^ i X, 
'o V 



This equation reflects the probabihstic parton cascade interpretation of the LLA |34, 35, p^ , 
in which the change of the gluon distribution on the left hand side is governed by the balance 
of gain (+) and loss (— ) terms on the right hand side. Notice that ( |38|) is free of infrared 
divergences, because the singularities in (^3|) at z = and z = 1 cancel between gain and 
loss terms. A diagrammatic illustration of these gain and loss terms is shown in Fig. 4. Also 
notice that the gluon fusion terms (second and fourth term) involve the 2-gluon distribution 
g^'^\xi, X2, k'j_), the presence of which causes the evolution equation to be non-linear. 

By means of a Mellin transformation the multiple convolutions of z-integrals inherent in 
the iteration of eq. (|3^ can be converted into products of independent successive interaction 
probabilities. Let me define the gluon distribution in the moment representation as 

g{u;,kl) := / dxx''g{x,kl) = / dye~^^ xg{x,kl) , (39) 

JO 70'--' 
where y = ln(l/x), in which the variable w is conjugate to ln(l/a;), implying that the low 
X behaviour is determined by small values of w. Analogously, the 2-gluon distribution is 
represented as 



g^^\LO,z,kl) :-- 



dxx^ g^'^'^ ( X, 

\ z 



-X, k\ 



(40) 



which carries an additional z-dependence. In general g^"^^ is a complicated correlation 
function not available in analytical form, except for certain special cases 37]. It is 
therefore inevitable on an analytical level to assume a phenomenological form of g^"^^ that 
allows to convert eq. (|38|) into a tractable linear form. This can be achieved with the 



following ansatz of product form |28 

1 



g^^^ I X, 



p{kl) g (x,k'i) g 



-X, k 



(41) 



where p{kj_) is a parameter that characterizes the magnitude of the probability for finding 
two gluons at the same point in phase-space, depending on their typical transverse size r± ~ 
l/k±. Since the 2-gluon correlation must become substantial when the number of gluons per 
unit area ng/{'KR\) becomes so large that the gluons overlap spatially {R_\_ 0.5 — 1 fm), 
one expects that p = 0(1) when k^ {Qq = 1 GeV), and monotonically increasing as 



14 



k"^ — > A^. Using (^) in (^), the 2-gluon distribution in the moment representation can 
be approximated in the soft hmit (z <^ 1) as 



g'^^\u;,z,kl) ^ p(kl)—g{u;,kl) 



(42) 



In the moment representation the evolution equation (^) now becomes an hnear alge- 
braic equation for the Mellin transformed gluon distribution: 



where 'y{uj,k'j_) plays the role of a generalized anomalous dimension, 



as(A;j 



2tt 

27r 



e{kl - Ql) 



A, 



g^ax 



(a;) 



kl 



(43) 



(44) 



The functions A{uj) are given by 



Ag^ggiuj) 



Agg^giuj) 



2Ca 



11 1 

12 ^ 



+ 



1 



1 



LO + l uj + 2 cj + 3 



5(a;) 



vr^ p 



A, 



LO 



A 



9-*gx 



A 



3 1 
2 



1 



vr^ p 
2uj 



uo + l uj + 2 
1 2 



2S{uj 
2 



(45) 

d^'n.T{zy\/ dz and 



_a; + l a; + 2 u; + 3 

where S{lo) = il){u) + 1) — V'(l) with the Digamma function il){z 
—'(/'(I) = "fE = 0.5772 the Euler constant. 

The anomalous dimension ^{uj,k\), eq. (|4^), reduces at /cj^ » to jqcd, the QCD 
anomalous dimension in the LLA ||3^, since in this kinematic region k{x) = 1 and therefore 
A;^ = 0. However, at k"^ ~ Qq, A^(/c^) becomes non-zero, so that the evolution of the 
gluon distribution receives modifications due to the coupling of gluons to the x-held. In the 
region Qq > kj_ > A'^, the perturbative QCD contributions oc Og vanish per construction, 
so that the gluons solely interact with the x-field, the coupling to which increases, because 
K ^ 0, i.e. A;^ 1. This behaviour is evident in Fig. 5, which shows 7(cj, k'j_) versus to for 
different values of /c^ . 
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The formal solution of eq. ( ^3[ ) is 

g{uj,kj_) = 



(46) 



from which the x distribution can be reconstructed by considering the inverse Mellin trans- 
form 

xg{x,kl) = ^ J^dux-'^g{uj,kl) , (47) 

where uj is now a complex variable and the contour of the integration C runs paralell to 
the imaginary axis. For the full anomalous dimension (^) this inversion must be done 
numerically [^]. 

3.2 Analytical estimates for a;-spectra and gluon multiplicity 

To exhibit the main features of the evolution of the gluon distribution in the presence 
of the x-field, it is instructive to make some analytic estimates. Of particular interest is 
the low x-region, because the soft, small x gluons are most preferably radiated, but at the 
same time take away only a very small portion of the total energy. For simplicity I will 
divide the kinematic range of the /c^-evolution into two distinct domains as indicated in 
Fig. 2: 

(i) Q"^ > k\> Qq and k^/K^ » 1: In this region A^^- ~ and dX^/dlnk^ ~ 0, so that 
(E3) reduces to 



7(0;, /Cj 



a. 



27r 



(48) 



(ii) Ql>k\> A2 and k\/K^ 1: Here A^^ ^ (47r)-i and d\^/d\nk\ > 0, in which 
case (113) becomes 



7(0;, /Cj 



27r 



A 



'A2 
k 



2 I ^^ga-^xi^) 



(49) 



As stated before, the low x region corresponds to the limit ^ 0, so that an expansion 
of 7(w,A:^) around uj = gives the dominant contributions at small x. Up to to order lo 
one has in the region (i) 



A: 



d 



\D.g{u),k'\ 



27r 



7«(u;,fci), 



(50) 
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whereas in the region (ii) one gets in the same order of approximation 



k 



d 



In g{uj,k] 



27T ^k^, 



»7r" 
3a; 



(51) 



The accuracy of the expressions 7*^*^ and 7^**-' in the relevant range of uj is exhibited in Fig. 
6, where the exact expression 7(0;, k'j_) is compared to the smah-w expansions ( |50D and ( |5l| ) 
at large and small k'j_. Evidently the approximation is rather good for lo < 2. 

Eqs. (H) and ® can now be solved analytically with the saddle point method. The 
/c5_-dependence of Os is given in (U) and for A-^^ = |^/(47r), eq. (|35|) , I will use here the 
form 

,2 ^ _ OiQl - kl) HQl/kl) 



(52) 



Att ln(Q2/A2) ' 

which has the required properties that = {^^^^ = 0) for k'j_ > Qq and A^^ (4vr)~-^ 
(.^j(. 1) for k"^ — > (as before denotes the Fourier transform of ^{x) in eq. 
Introducing the variables t for the region (i) and u for the region (ii), 



t{kl) 
u{kl) 



?2 J, '2 



dk^ as{k^ 



2tt 



2TTb 



In 



4dk^ \{kl) A2 



ln(QVA2) ' 
ln(A:2/A2) 

p 1 - {Ql/kD ( 1 - HQi/kD) 



Ikl k'^ \ 27r '^klj 87r2 {Ql/ K^)\ii{Ql/ k-^) 

the combination of (|4^ ) and (^) yields for the kinematic domains (i), respectively (ii): 



(53) 



xg{x,t) 



: / da; exp[c^y + i/«(a;)tl 5(0., 0) , i/«(a;) = 2C7^ f- " ^"l 

27n Jc L J ^tj 12/ 

j;ff(x,ii) = — / da; exp + i/(**)(a;)M ^(w, 0) , i/(**)(w) = — , (54) 
ZiTi Jc" L J 3u; 

where y = ln(l/x). The saddlepoint a;^*'' of the integrand in ( ^^ is determined by the 
condition 

duj 1^ J ,„=,„(') 







(55) 



and similarly o;^ . Using the method of steepest descent then gives the following results. 



(i) In the region Q > k^ > Qq. 



xg{x,k\) 



Q^>kl>Ql 



Ni{x) + Gi{x,k'i) exp 



ll^A.. „2 



127r6 



Hiiki) 



exp 



7rb 



Hi{ki)\n 



Ni{x) = xg{x,kl=Q^) = 6(1 - x)6ikl - Q'^) 
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In 



ln(QVA2) 



ln(A:2/A2) 
(ii) In the region Qq > kj_> Q^: 



xg{x,k\) 



1/4 



ln( i 

X 



-3/4 



Q2>fc2 >A2 



N2{x) - G2(x,fci) exp 



'f«.(.i).n(i 



iV2(x) 

G2(x,fei) 
H2{kl) 



xg{x, k\ = Ql) 



^H.ikl 



1/4 



ln( - 

X 



-3/4 



(57) 



1 - {Ql/kl){l-HQl/kl)) 



{Ql/A^)HQl/A^) 

In Fig. 7 the x-spectra of ( p6| ) and ( |57[ ) are shown for different values of k'j_ with fixed 
Qo = 1 GeV. Two different initial distributions were chosen to start the evolution from 
k"^ = (3.5 GeV)^, one flat in rapidity and the other one a Gaussian form. The parameter 
p introduced in ( pi] ) was set equal to one. For k'j_ > Qq the gluon distribution xg{x, k'j_) is 
just the well known LLA solution with its strong increase at small x as kj_ decreases. For 
k'j_ < Qq however, the effect is reversed such that g{x,k'j_) decreases as k'j_ falls below Qq. 
This suppression, which is particularly substantial at small x, reflects the "condensation" 
of gluons in the collective background field x- 

The /c^-dependence of the total gluon multiplicity can also be estimated in the above 
approximation. The gluon multiplicity is given by the lo = moment, 

»i 



Ng{ki) = g{co = 0,ki) 



dx g{x, kj_) . 



(58) 



Using eqs. (^3|)-(p^), one arrives in the soft limit {z <^ 1) |3^ at the following integral 
equations that govern the approximate behaviour of the gluon multiplicity for the two cases 
(i) and (ii): 



ki 



{ii) 



d 



d 



Ngikj 



.^N.ikl) 



^0 dk 



N,{k, 



(59) 
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On account of the /c5."dependence of Og ( p4[ ) and A^^, (52), the corresponding solutions are 
obtained as: 



(i) Ng{k 



N, 



, /ln(QVA2)\-V4 exp [2 VICaM) MQ^A^) 
W j i„/7 9 /A9\ r r : — : — — r- 



iii) N,{ki) ^^^^^^^^ = NM) exp 



pk^ l + [l + Hkl/Ql)f - 2kl/Ql 
2 kl ln(Qg/A2) 



(60) 



In the region (i) > /c^ > Qq; multipUcity coincides with the QCD result |33], 
characterized by a rapid growth as the gap between the hard scale and k\ increases. 
On the other hand, in the region (ii) > k\ > h?, the multiplicity becomes strongly 
damped. The exponent is always negative so that the number of gluons rapidly decreases 
and vanishes at k\ = 0, ensuring that no gluons and therefore no color fluctuations exist 
at distances R>^Ar^. This behaviour is evident in Fig. 8, where Ng{k\) is plotted versus 
/k\, starting from A^/Q^ < 1. 

It must be emphasized that the perturbative evolution of gluons is cut off at Qq, and in 
the transition region Qq> kj_ > A^ the evolution is purely non-perturbative, although it is 
described here as an extension of the perturbative evolution above Qq and treated on the 
same footing. 

3.3 Flux tube configurations of gluons interacting with the mean field x 

In Sees. 3.1 and 3.2 the evolution of the gluon configuration between the fragmenting qq 
pair was analyzed in terms of the non-perturbative modification as a function of ^xi^i.) ~ 
^^/(47r). However, the coupling strength or between the gluon field and the x- 
field must be determined by the dynamics itself, since is the Fourier transform of the 
X-dependent coupling function £,{x) in eq. (^). 

The dynamics is governed by the the coupled system of equations eqs. ([T9|)-(^l|), which 



can now be solved numerically by utilizing the definition (36) together with the general 
solution for the gluon distribution given by ( p6| ) and (p7[). To do so, one obtains the 
expectation value of eq. (|20|) by multiplying with the multigluon state vector (P| and 
\P) from the left and right, respectively. In the mean field approximation, x = X is a c- 
number function, so that this operation affects only the F^^F^^ term which gives the gluon 



distribution by virtue of ([Sq) . With the explicit solution of the gluon distribution (|5q) and 
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( p7D inserted, the remaining task is to solve the single equation (^) for the x-field. (Recall 
that the equation (^) for the [/-field couples only to x ^'^'^ is in principle readily solved 
once the solution for x is known.) 

An interesting phenomenological application 22, 41 1 is to calculate the string tension 
t, which characterizes the linearly rising potential between the qq-pair due to the gluon in- 
teractions, usually obtained by fitting heavy quarkonium spectra fi^ with a non-relativistic 
potential of the form 

VUr) = + tr\ (61) 

where a = 40^/3. Typical fit values for the string constant t range from 750 MeV/fm to 
950 MeV/fm. 

Here I shall estimate the string constant within the present approach on the basis of the 
equations of motion ([T9|)-(pl|). Similar calculations have been done earlier in the framework 
of the static MIT bag model |]4^ and the Friedberg-Lee soliton model 41]. 

In the following I will consider the fragmentation of a heavy qq pair within an adiabatic 
approximation. That is, a quasistatic treatment is employed which neglects the motion 
of the qq pair and considers the instantanous gluon configuration in between the pair. 
This should provide is a reasonable approximation, because one can view the z*'* gluon as 
being emitted from the qq pair plus gluons 91,92,- ■■ ,gi-i, with the spatial coordinates 



of these "sources" being frozen during the emission of the gluon i [44|. In a space-time 
picture of the fragmentation of the qq and its emitted gluons it is the change of the typical 
transverse momenta k± or transverse separation r± oc kj_^ of gluons which governs the 
dynamical transition from short distance, unconfined regime to long distance, confined 
stage because there must be a critical separation of color charges beyond which the 
total color is screened. In the present case, the role of this non-perturbative phenomenon is 
played by the x-held. In the adiabatic approximation, there is no explicit dependence on the 
longitudinal variables here, because the separation in the transverse plane is independent 
of when during the cascade, or where along the jet axis, a gluon was produced p5| . 

Thus, within the adiabatic treatment, one may use the separation Rqq of the receding 
qq as a measure for the typical gluon transverse momenta ~ Rq^- By minimizing the 
energy per unit length of this system, one obtains then for each gluon configuration at a 
given Rqq the form of x a-nd the string tension t. The total energy per unit length, the 
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string tension t, receives contributions from both the x-held in the low energy regime and 
the gluon field in the high energy domain, which on account of the the trace anomaly (^) 
is given by: 



E 



qq 



dA 



1 



+ 



J dA k{x) {P 



F,u,aFr\P) . (62) 



Here A is the cross-sectional area of the flux tube of the gluons between the qq, perpendicular 
to Rqq. In one loop order the beta-function is (3 = —ba'^ with b = 33/(127r) for Nj = 0. 
Then, by using eq. ( |36| ) to express the second integral in terms of the gluon distribution 
g{x, k'j_), assuming cylindrical symmetry along the Rgg axis, and minimizing t with respect 
to X, one arrives at the following non-linear integro-differential equation (r is the radial 
coordinate perpendicular to Rqq)' 

dVix) b 



d2 1 d 

dr"^ r dr 



X{r) + 



dx 



+ {P+ f-Ig{Rqq) J drrnix) = 0, 



(63) 



where 



hi^m) = / d^k^as{kl) / dxxg{x,kl). (64) 

Jl/Rqq Jo 

In pulling Ig out of the r-integral, it is assumed that the spatial distribution of gluons 
is approximately homogenous. Eq. ( |63[ ) determines the energetically most favorable flux 
tube configuration. However, a physical meaningful flux tube solution has to satisfy the 
constraint that the system as a whole, qq plus gluons, must form a global color singlet, 
implying that all of the color flux that originates from the q must be directed towards the 
q. In other words, the total color electric flux through a plane between the q and q must 
be equal the color charge Qq = QsTa on one of them |22, 41]. This translates into the 
requirement that the gluons in the flux tube streched between q and q with certain Rqq 



carries a total color charge squared that is equal to Q'i, the one of qq. Define 



A 



j d\ Jg{Rqq) , 



(65) 



dx g{x, kj 



(66) 



where A is the cross-sectional area of the flux tube, and 

Jg{Rqq) = / d^ki 

Jl/Rqq JO 

is the total number of gluons radiated from the initial point of (/^-production up to Rqq. 
Then the above constraint then reads 

levr 



^9 = Ql = idTa-Ta] 



as{k 



k^=RZ 



(67) 
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Combining eqs. (pq)-(|67D, one arrives at 
' Id 



dr"^ ^ r dr 



X{r) + ^ + fbasiR-?) / drr.ix) = , (68) 



which is now independent of the overall boost momentum (P^), and thus of the initial hard 
scale = (P^)^, as it should be. 

Solving eq. ( |68D numerically I^Bj, subject to the the boundary conditions x'(0) = 
and x(oo) = Xo yields the solutions for t, x ™d shown in Fig. 9. The reasonable 



parameter values |17| xo = f-w and bag constant B = (150 MeV)^ were chosen for V{x)^ 
eq. (0), and the coupling function k was taken of the form ([ll[). One sees that with 
increasing separation Rqq of q and g, the gluons first multiply which results in a growing 
string tension, but then gluon condensation sets in, yielding a saturating behaviour of t 
with the string constant approaching f ~ 1 GeV/fm (Fig. 9a). For Rqq ~ 1 fm the 
gluon field is completely confined within a fiux tube of radius r ~ 1 fm (Fig. 9b). For 



comparison, a simple estimate within the MIT bag model [43| gives t = 910 MeV/fm, but 
a rather large tube radius of 1.6 fm. Detailed calculations with in the soliton model p^ ] 
gave qualitatively similar results. Finally, Fig. 10 shows in correspondance to Fig 9b the 



form of the potential V{x) = V{x, C/)|mq=0) defined in eq. (14), and the effective squared 
mass M^(x) = d^V {x) / dx^ for the same parameter values as above. 

4. SUMMARY AND OUTLOOK 

The effective QCD field theory approach presented here to describe the dynamics of high 
energy partons in the presence of a collective confinement field provides a framework that 
has the potential to be developed towards a systematic description of the hadronization 
mechanism. The corresponding effective action has been constructed such that it 

(i) incorporates both parton and hadron degrees of freedom; 

(ii) recovers the exact QCD (Yang-Mills) action with its symmetry properties at short 
space-time distances; 

(iii) merges into an effective low energy description of hadronic degrees of freedom at large 
distances; 
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(iv) allows for a dynamical description of parton-hadron conversion on the basis of the 
resulting equations of motion. 

As an exemplary demonstration, the approach was applied to the evolution of a fragmenting 
qq pair with its generated gluon distribution, starting from a large hard scale all the 
way downwards to A^. The transformation of the initially high virtual gluons to a gluon 
condensate field x was studied in terms of the coupled evolution of the gluon distribution 
and the mean field x- The solution of the equations of motion yields color fiux tube 
configurations with an associated energy per unit area (string tension) of about 1 GeV/fm, 
consistent with the common estimates. 

In perspective, important points to be addressed in the future, are: 

(i) The establishment of the relation with the exact renormalization group equation for 
the effective action as derived by Renter and Wetterich is desirable. This would 
allow to quantify the effect of consecutively integrating out all quantum fiuctuations 
of gluons and quarks with momenta larger than some infra-red cut-off scale Qq, the 
variation of which determines the confinement dynamics. 

(ii) With the inclusion of quark degrees of freedom and possibly quantum fiuctuations of 
the X- [/-fields, one could calculate e.g. the mass spectrum of glueball and meson 
excitations as physical hadrons. This would provide a complete description from a 
physical initial state, via a not directly observable deconfined partonic stage, up to 
the formation of observable hadronic excitations. 

(iii) Ultimately one would like to address the microscopic dynamics in full 6-dimensional 
phase-space with explicit inclusion of the color degrees of freedom and the local 
color stucture. This could be realized in a transport theoretical formulation similar 



as in Ref. |4S], in which the partons propagate with a modified propagator that 
embodies the effects of the mean field x in the effective mass. As the confining field 
becomes significant the effective mass increases and asymptotically becomes infinite 
so that the propagation of color fiuctuations ceases. 

(iv) The possible applications are manifold. One particular interest is the expected 
(non) equilibrium QCD phase transition in high energy systems as in heavy ion col- 
lisions or the early universe, an issue which could be addressed along the lines of 
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Campbell, Ellis and Olive [17| in combination with the space-time evolution of the 
multi-parton system |47| in the presence of the collective field. 
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APPENDIX A 



For completeness a "derivation" of the evolution equation (pq) is given in the spirit of 



Lipatov |49|. The full propagator of a single gluon of momentum k may be represented as 



A;2(l + n(fc2)) l + n(A;2) 



(69) 



where Z^JJ^ = di^ty/k"^ is the free propagator, and with the choice of gauge for the gluon 



fields 7? • ^ = 0, 



kfj,Tji, + kui]^ 
k ■ r] 



df^f^ik) 



(70) 



and k^^d^u{k)\k'i=Q = guarantees that only 2 physical (transverse) gluon polarizations 
propagate on mass shell. The self-energy part 



n(fc2 



n 



9^9 



{k^) + Tlg^^{k 



(71) 



contains both the gluon self- interactions and the "medium" corrections due to the coupling 
to the confining background field x Expanding H^^g and H^^^ in powers of the squared 
couplings g1 and respectively, the contribution to one loop order is determined by the 
total gluon "decay" probability, i.e. the probability of losing a gluon out of a momentum 
space element between A;^ and k"^ + dk'^: 

d 



Wg{k 



ain(A;'2/A2) 



n 



'J 2 



9^9 \ 



) + n,^^(/c 



/2^ 



9^9 



+ W 



9^X ■ 



(72) 



Here, Wg^g and Wg^^ are the inclusive probabilities for a gluon to emit (absorb) another 
gluon, due to the self-interaction, respectively the interaction with the x-held, corresponding 
to the diagrams in Fig. 4, 



Wg^g{k^) 
Wg^xik^) 



dx 



dx 

Jx X' 



1 dx' r 

x' 
1 dx' 

./ 



^9^99^"^ , X, A; ) + Wgg^g(x ,x,k ) 
^g—*9xi'^ , a;. A; ) + Wgg—^^{x ,x,k ) 



(73) 



The individual contributions in square brackets w{xi, X2, k'^) can be obtained in the stan- 



dard fashion [pQ, 51| by evaluating the cross-section ratios (c.f. Appendix B) 



/fc2 daW 



a 



(0) dzdk'^, 



87r2 



(74) 
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where g denotes the appropriate couphng for the process under consideration (here g = gs 
or g = ^^), o"^'^) is the Born cross-section for the production of a gluon a and a^^^ represents 
the first order correction associated with the "decay" a be. 



For the process g ^ gg and its reversal gg ^ g the probabihty distributions (33) are 
well known |]50| , |38| . Assigning the momentum fractions as xi — > X2, {xi — X2) for g ^ gg 
and xi, (x2 — xi) — > X2 for gg g, one has 



Wg^gg{Xl,X2,k 



Wgg^g{Xl,X2, k^) 



2-K 



X2 



2tt 



h? Xi{x2-Xl) 



^9-^99 



(75) 



where in the second expression the factors in square brackets arise from the difference 
of phase-space and flux factors for fusions compared to branchings. The color factor is 



"^99^9 



1/8 and 



-fg^ggiz) = 2Ca [z{1-z) + 



+ 



1 



(76) 



1 — z z 

where Ca = -^c = 3, and z is the fraction of x-values of daughter to mother gluons. 

The new, additional processes are the friction process g gx, corresponding to energy- 
momentum transfer from gluons to the x-field, and the fusion process gg by which 
two gluons couple color neutral to the x-field and "annihilate" 0. As outlined in Appendix 
B, one arrives at (with the assignment xi — > X2, (xi —X2) for g ^ gx and xi, (x2 — xi) X2 
for gg x) 



w. 



where c, 



9— »ffx('^i' -^Si k ) 
1/8, and 



■9g^x 



Mil 
27r 

_ 

2-K 



'~19^9X 



57r C. 



99-*9 



1 ( 1 + Z^ ^ 
4 i 1 - z 



A^ Xl(x2-Xi) 



^X^99 



(77) 



z + 



1 



(78) 



^Since the conversion of partons into hadrons in the process of fragmentation is an irreversible process, 
the spontanous production of gluons by the x-field x ^ 99 y ^ well as the energy transfer from the x-field to 
the gluons, gx ~* g, are omitted. These latter interactions would counteract the transition, which certainly 
is possible in the sense of local fluctuations, but globally, and in the average, the parton-hadron conversion 
is a one-way process in the present context. 
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The total interaction probability Wg (72) determines via the unitarity condition 
the Sudakov formfactor Fg, 



Fg{Q^,k^) = exp 



fc2 k 



— Wg{h!'^) 



(79) 



which is the probability that a gluon does not at all interact (i.e. emit or absorb other 
gluons) while degrading its virtuality from to k"^. 

The self-energy part (71) is now readily evaluated on the basis of eq. ([7^), and inserted 
in the representation (|69|), one obtains the single gluon propagator at one loop order, 



D^.u{k) = Dl{k) 



^ + / ~rT / —^{^ ^x,k 
Jk^ k ^ Jx x' 



(80) 



The corresponding "jet calculus" |27] Greens function L'g( Q'^) of eq. (p^) describes 

how a system of gluons evolves in the variable x and the virtuality k"^ through the gluon 
self- interactions and in the presence of the confining background field x- It is given by 
the convolution of the single gluon propagator ( |69[) with the gluon distribution function 
g{x,k^). Defining 

Dg{x,k'';xo,Q^) = Qf^uD'^-'ik) gix,k'') , (81) 

the self-consistent iteration of one loop contributions to all orders within the LLA gives the 
evolution equation for the gluon distribution with respect to the variables x and k^: 



^'^'O^ 9{x,k^) = k^-g^ Dg{x,k'^;xo,Q'^^ 



+ 



1 dx' 

X 

A2 



+ 



2tt 



\{k'^) r /•! dx' 



\g{x\k^)lg^ggl^^}j " ]^g{x,k^) dz jg^gg (z) 



k 



fc2 



dx' g^'^^ (x, x', k^) Tgg^g (x, x', X + x') 



1 

- / (ix'5(^)(x - x',x',fc^) r,,,,^^ (x - x',x',x)) 



fc2>Q2 



2-K 
27r 



r' / X \ 

yg{x',k'^)-fg^g^\-;j - g{x,k^) dz-fg^g^ 
8^^^ ) [J^ ^^'5^^^ k^) ^gg^x i^^ 2;', X + x') I 



where the factor 1/2 in the first (third) term arises from the indistinguishability of the two 
gluons emerging from (coming in) the branching (fusion) vertex. The function g^'^^ (xi, X2, fc^) 
denotes the 2-gluon density, and the gluon fusion functions T are defined in accord with 
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and (|73) as 

i 12^3(a;i,X2,X3j = Ci2^3 — 73^12 [ I = Cl2->3 2~ ">'3^21 I • (83) 

Xn \X3/ Xn \X^J 



^3 -^3 



with X3 = xi + X2. Changing to variables {x,xi,k ) {x,z,kj_) and using p^), one 
immediately arrives at eq. (|3l|). 

APPENDIX B 

Here I will outline the explicit calculation of the interaction probability densities 'Jg^gx 
and Tgg^^ that appear in the evolution equation ( ^2|) or ( ^8|) in addition to the usual 
probability densities 'Jg^gg and Tgg^g. Notice that the vertex corresponding to three gluons 
coupling to the x-field is unphysical and therefore to be excluded, because x is required to 
be a color singlet field. On the other hand, the coupling of four gluons to x is possible, 
however in the LLA such diagrams are kinematically suppressed and can be neglected Q. 

Let 0"^^^ denote the spin and color averaged cross-section for the production of a gluon 
at order in perturbation theory. The probability distribution 'ja^bc in the variable 
z = Xb/xa for the emission of a gluon b in the process a — > be, is the given by the ratio of 
cross-sections 

^ ^ = 8^ ^^-'^^'^ ^ ' 
where g is the appropriate coupling of the process, cr^'^^ is the lowest order cross-section for 
the production of a gluon a and a^^^ represents the first order correction associated with 
the "decay" a be. The vertex function associated with the general ggx coupling is easily 
obtained from the interaction Lagrangian C[ip,A,x], eq. (0), as: 



ab 



{ki ■ k2)g^u - (1 - a) kif,k2u 



(85) 



where Cx(^^) denotes the Fourier transform of ^(x) in coordinate space, ki, k2 are the gluon 
momenta and the convention is that all four-momenta are directed into the vertex. The 
process g ^ gx gives then (setting a = 1): 

where 



16 , , ,, 

a,a',b,b' si,S2 
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where the factor 1/16 in front results from the averaging over initial 2 transverse polariza- 
tions and 8 color degrees, and it is summed over final color and spin polarizations Sj. The 
sum over gluon polarizations si,S2 must be performed over transverse polarizations only. 
This is achieved by the projection 



2 ■ 



Assigning the momenta = {k^,k ,k±) of incoming (outgoing) gluon ki (^2) and the 
momentum k transferred to the x-^eld as 



k 



1 



k+ 




(1- 


-z)kt, 


zkf 


k^ 

' zkf 



' = . (89) 



and carrying out the appropriate change of integration variables, the result is 

r''i dk 



1 / z 
- l-z + 2- 



4 V l-z 



(90) 



Hence, one can read off 



cia« 1 fl + Z^^ 



In complete analogous manner the process gg ^ X can be calculated. The procedure 
is to evaluate x ~^ gg^ with incoming momentum k and the outgoing momenta ki and 
k2. Using the formula ( |83| ) one obtains the 2-gluon fusion function for the reverse process 
gg — > X- The result is: 

Tgg^^iXl,X2,X3) = Cgg^^^^-f^^gg(^) , (X3=Xl+X2), (92) 
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FIGURE CAPTIONS 



Figure 1: Typical shape of V{x,U), eq. ([T^), with dV/dx = at x = Xo and V = B at 
X = 0, where B = 6xo/4 is the bag constant. 

Figure 2: Schematics of the parton shower evolution of a fragmenting qq pair with its 
gluon configuration as the virtualities of the partons gradually degrade, starting from the 
hard scale Q^. At large gluon virtualities /c^ the shower develops by perturbative branching 
processes, but at k"^ ~ Qq non-perturbative fusion and friction processes set in, such that 
at k'^ = no colored fluctuations remain. 

Figure 3: Diagrammatic representation of the two-point Greens function of gluons, in- 
cluding both the gluon (self) interactions and the effective interaction with the confining 
background field x (indicated by the dashed lines). This gluon propagator describes the 
evolution of a gluon from a chosen cascade branch in x and k^, starting from xq and Q^. 

Figure 4: Diagrams for the interaction probabilities that determine the evolution of the 



gluon distribution according to (38). In the probabilistic interpretation terms with positive 
signs leadt to a gain of gluons and terms with negative sign to a loss. The gluon with 
momentum fraction x is the "observed" particle. 

Figure 5: The anomalous dimension 'y{uj,k'j_) of eq. (^) versus lu for different values of 
kl. 

Figure 6: Comparison between exact expression for the anomalous dimension (p^ ) and 
the approximations ( |50| ) [top] and (p|) [bottom]. 



Figure 7: Behaviour of the distributions xg{x, k'j_) of ( |56D and ( [57| ) for different values of 
kj_. Top part shows result for a flat initial distribution xg{x, Q"^) = a(l — xY with a = 2.8, 
c = 5.3, and bottom part for a Gaussian initial distribution xg{x, Q^) = 1/V 27rc^ exp [— (x — 
d)/{2c^)] with c = 0.3, d = 0.5. (Q = 3.5 GeV, Qq = 1 GeV, A = 0.23 GeV). 

Figure 8: Evolution of the gluon multiplicity Ng{k'^) from down to A^. At first the 
gluons multiply, but at /c^ ~ Qq a condensation sets which is complete at A^ [Qq = 1 GeV, 
A = 0.23 GeV). 
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Figure 9: a) String tension t versus separation Rqq of the qq pair, and b) the solutions 
for X f-iid K at Rqq = 1 fm versus r which is the radial coordinate perpendicular to Rgg 
(Xo = U and B^^ = 150 MeV). 

Figure 10: a) Form of the potential l^[x(?')] and b) the effective squared mass M^[x(r)] 
of the X field, both at Rqq = 1 fm. (xo = fn and = 150 MeV). 
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